Vortex state in a doped Mott insulator 
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We analyze the recent vortex core spectroscopy data on cuprate superconductors and discuss what 
can be learned from them about the nature of the ground state in these compounds. We argue 
that the data are inconsistent with the assumption of a simple metallic ground state and exhibit 
characteristics of a doped Mott insulator. A theory of the vortex core in such a doped Mott insulator 
is developed based on the U(l) gauge field slave boson model. In the limit of vanishing gauge field 
stiffness such theory predicts two types of singly quantized vortices: an insulating "holon" vortex 
in the underdoped and metallic "spinon" vortex in the overdoped region of the phase diagram. 
We argue that the holon vortex exhibits a pseudogap excitation spectrum in its core qualitatively 
consistent with the existing experimental data on Bi2Sr2CaCu208- As a test of this theory we 
propose that spinon vortex with metallic core might be observed in the heavily overdoped samples. 



I. INTRODUCTION 

Nature of the ground state as a function of doping 
remains one of the recurring unresolved issues in the the- 
ory of high-T c cuprate superconductors. The problem is 
partly due to formidable difficulties related to the the- 
oretical description of doped Mott insulators and partly 
due to experimental hurdles in accessing the normal state 
properties in the T — > limit because of the interven- 
ing superconducting order. Probes that suppress super- 
conductivity and reveal the properties of the underlying 
ground state are therefore of considerable value. So far 
only pulsed magnetic fields ffl in excess of H c2 and im- 
purity doping beyond the critical concentration B have 
been used towards this goal. Here we argue that the 
vortex core spectroscopy performed using scanning tun- 
neling microscope (STM) can provide new insights into 
the nature of the ground state in cuprates. We analyze 
the existing experimental data and conclude that 

they imply strongly correlated "normal" ground state, 
presumably derivable from a doped Mott insulator. We 
then develop a theoretical framework for the problem of 
tunneling in the vortex state of such a doped Mott insu- 
lator. 

In the vortex core the superconducting order parame- 
ter is locally suppressed to zero and the region within a 
coherence length £ from its center can be to the first ap- 
proximation thought of as normal. Spectroscopy of the 
vortex core therefore provides information on the nor- 
mal state electronic excitation spectrum in the T — ► 
limit. More accurately, the core spectroscopy reflects the 
spectrum in the spatially non-uniform situation where 
the order parameter amplitude rapidly varies in response 
to the singularity in the phase imposed by the external 
magnetic field. In order to extract useful information 
regarding the underlying ground state from such mea- 
surements a detailed understanding of the vortex core 
physics is necessary. So far the problem has been ad- 
dressed using the weak coupling approach based on the 
Bogoliubov-de Gennes theory generalized to the <i-wave 



symmetry of the order parameter and semiclas- 

sical calculations The early theoretical debate 

focused on the existence or absence of the vortex core 
bound states 0|u||n|. This debate, now resolved in fa- 
vor of absence of any bound states in pure d x - y i state 
P, ^0|p6| , has somewhat eclipsed the possibly more im- 
portant issues related to the nature of the ground state 
in cuprates. 

The body of work based on mean field, weak cou- 
pling calculations l8|- |l0Ul^] yields results for the local 
density of states in the vortex core which exhibit two 
generic features: (i) the coherence peaks (occurring at 
E = ±Ao in the bulk) are suppressed, with the spectral 
weight transferred to a (ii) broad featureless peak cen- 
tered around the zero energy. Here we wish to empha- 
size the heretofore little appreciated fact that these fea- 
tures are qualitatively inconsistent with the existing ex- 
perimental data on cuprate superconductors. STM spec- 
troscopy on Bi 2 Sr 2 CaCu208 (BSCCO) at 4.2K indicates 
a "pseudogap" spectrum in the vortex core with the spec- 
tral weight from the coherence peaks at ±Ao ~ 40meV 
transferred to high energies, and no peak whatsoever 
around E = Q B. Recent high resolution data on the 
same compound |(| confirmed these findings down to 
200mK and found evidence for weak bound states at 
±7meV. Experiments on YBa 2 Cu 3 7 (YBCO) [§ also 
indicate low energy bound states, but are somewhat more 
difficult to interpret because of the high zero-bias conduc- 
tance of unknown origin appearing even in the absence 
of magnetic field. 

The fundamental discrepancy between the theoretical 
predictions and the experimental findings strongly sug- 
gests that models based on a simple weak coupling theory 
break down in the vortex core. The pseudogap observed 
in the core hints that the underlying ground state re- 
vealed by local suppression of the superconducting order 
parameter is a doped Mott insulator and not a conven- 
tional metal. Taking into account the effects of strong 
correlations appears to be necessary to consistently de- 
scribe the physics of the vortex core. Conversely, study- 
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ing the vortex core physics could provide information es- 
sential for understanding the nature of the underlying 
ground state in cuprates. 

The first step in this direction was taken by Arovas et 
al. ]rj| who proposed that within the framework of the 
SO (5) theory vortex cores could become antiferro- 
magnetic (AF). They found that such AF cores can be 
stabilized at low T but only in the close vicinity of the 
bulk AF phase. In contrast, experimentally the pseudo- 
gap in the core is found to persist into the overdoped 
region jij . More recently microscopic calculations within 
the same model jl9| revealed electronic excitations in 
such AF cores with behavior roughly resembling the ex- 
perimental data. Quantitatively, however, these spectra 
exhibit asymmetric shifts in the coherence peaks (related 
to the fact that spin gap in the AF core is no longer tied to 
the Fermi level) not observed experimentally. These dis- 
crepancies suggest that generically cores will not exhibit 
the true AF order. Finally, these previous approaches 
are still of the Hartree-Fock-Bogoliubov type and can- 
not be expected to properly capture the effects of strong 
correlations. 

Here we consider a model for the vortex core based 
on a version of the U(l) gauge field slave boson theory 
formulated recently by Lee |2(|. Originally proposed by 
Anderson the slave boson theory was formulated to 
describe strongly correlated electrons in the Cu02 planes 
of the high-T c cuprates. Various versions of this theory 
have been extensively discussed in the literature p2|-p6| . 
Interest in spin-charge separated systems revived recently 
p9|,p8 29 2(J due to the realization that it provides a nat- 
ural description of the pseudogap phenomenon observed 
in the underdoped cuprates. The common ingredient in 
these theories is "splintering" of the electron into quasi- 
particles carrying its spin and charge degrees of freedom. 
Within the theories based on Hubbard and t-J models 
this splintering is formally implemented by the decom- 
position of the electron creation operator 



„t 



(1) 



into a fermionic spinon and bosonic holon bi. The 
local constraint of the single occupancy b\bi + f} a fia = 1 
is enforced by a fluctuating U(l) gauge field a. The 
mean field phase diagram is known to contain four phases 
distinguished by the formation of spinon pairs, = 
(eerer' fi a fja')) an d Bose-Einstein condensation of the in- 
dividual holons b — (bi) p6| , and is illustrated in Figure 

The effects of magnetic field on such spin-charge sepa- 
rated system is most conveniently studied in the frame- 
work of an effective Ginzburg-Landau (GL) theory for 
the condensate fields A and b. The corresponding ef- 
fective action can be constructed [BO 31 based on the 



requirements of local gauge invariance with respect to 
the physical electromagnetic vector potential A and the 
internal gauge field a: 




FIG. 1. Schematic phase diagram of the system with 
spin-charge separation in the doping— temperature plane, as 
applied to cuprate superconductors. 



/gl = |(V - 2za)A| 2 + r a | A | 2 + -u A \A\ 4 

+ |(V - ia - ieA)b\ 2 + r b \b\ 2 + ^u b \b\ 4 + v\A\ 2 \b\ 2 



— (Vx A) 2 + / gaugc . 
on 



(2) 



The factor of 2 in the spinon gradient term reflects the 
fact that pairs of spinons were assumed to condense, 
/gauge describes the dynamics of the internal gauge field 
a. We note that unlike the physical electromagnetic field 
A the gauge field a has no independent dynamics in the 
underlying microscopic model since it serves only to en- 
force a constraint. Sachdev [Q and Nagaosa and Lee 
fjlf assumed that upon integrating out the microscopic 
degrees of freedom a term 



/gauge — —(V x a) 



(3) 



is generated in the free energy. They then analyzed vor- 
tex solutions of the free energy (g) and came to the 
conclusion that two types of vortices are permissible: a 
"holon vortex" with the singularity in the b field and a 
"spinon vortex" with the singularity in the A field. Be- 
cause holons carry electric charge e the holon vortex is 
threaded by electronic flux quantum hc/e, i.e. twice the 
conventional superconducting flux quantum $ = hc/2e. 
Spinons on the other hand condense in pairs, and the 
spinon vortex therefore carries flux $o- Stability anal- 
ysis then implies that spinon vortex will be stable over 
the most of the superconducting phase diagram, while 
the hc/e holon vortex can be stabilized only in the close 
vicinity of the phase boundary on the underdoped side 
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©0. This is a direct consequence of the fact that 
singly quantized vortices are always energetically favor- 
able nil. 

As far as the electronic excitations are concerned, the 
spinon vortex is virtually indistinguishable from the vor- 
tex in a conventional weak coupling mean field theory: 
the spin gap A, which gives rise to the gap in the elec- 
tron spectrum, vanishes in the core. Consequently, the 
vortex state based on the results of Sachdev-Nagaosa- 
Lee (SNL) theory |3^^l[] does not exhibit the pseudo- 
gap in the core and suffers from the same discrepancy 
with the experimental data as the weak coupling theories 
|t| |To| based on the conventional Fermi liquid description. 
Moreover, no evidence exists at present for stable dou- 
bly quantized holon vortices predicted by SNL. What is 
needed to account for the experimental data is a singly 
quantized holon vortex stable over the large portion of 
the superconducting phase in the phase diagram of Fig- 
ure |l|. In the core of such a holon vortex the spin gap A 
remains finite and leads naturally to the pseudogap ex- 
citation spectrum. In what follows we show that under 
certain conditions the free energy (|^) permits precisely 
such solution. 

The results of the SNL theory are predicated upon the 
assumption that the "stiffness" a of the gauge field is 
relatively large and that singular configurations in which 

V x a contains a full flux quantum through an elemen- 
tary plaquette are prohibited. Consider now a precisely 
opposite physical situation, allowing unconstrained fluc- 
tuations in a. This amounts to the assumption that the 
/gauge term (||) can be neglected in (||), i.e. a — > 0. Phys- 
ically this corresponds to the "extreme type-I" limit of 
the GL "superconductor" (^) with respect to fluctuations 
in a. Based on Elitzur's theorem Jm) Nayak recently 
argued that the exact local U(l) symmetry of the model 
cannot be broken, implying absence of the phase stiffness 
term (^) at all energy scales. Our assumption therefore 
appears reasonable and in Section III. we shall give a 
more thorough discussion of the significance of the /gauge 
term for the vortex solutions of interest here. For the 
time being we shall assume that /gauge can be neglected 
and explore physical consequences of the resulting theory. 

/gl given by Eq. (||) is quadratic in a and with the 

V x a term absent the gauge fluctuations can be trivially 
integrated out. Within the closely related microscopic 
model this procedure has been recently implemented by 
Lee [EOl . The resulting effective free energy density reads 



/ /amp 
1 



2 2 
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^(V(f>-2V9 + 2eAY 
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(V x A)', 



(4) 



where we have set A = p&e 1 ^, b = p^e 1 ® , and 
/amp = (Vpa) 2 + r A p 2 A + ^uaPa 

+ (Vpf,) 2 + npl + -u b pf + v P \p\ 



is the amplitude piece. The most important feature of the 
effective free energy (||) is that it no longer depends on 
the individual phases <f> and 9 but only on their particular 
combination 



26>. 



(6) 



(5) 



Since the physical superconducting order parameter \E' = 
A*b 2 — pi\p\e~ lly ^~ 2d ^ it is reasonable to identify Q with 
the phase of a Cooper pair. Physically, the unconstrained 
fluctuations of the gauge field in Eq. (Q) resulted in par- 
tial restoration of the original electronic degrees of free- 
dom in Eq. (0). In the underlying microscopic model 
this means that on long length scales spinons and holons 
are always confined, in agreement with Elitzur's theorem 
Q,0. On lengthscales shorter than the confinement 
length, such as inside the vortex core, spinons and holons 
can still appear locally decoupled. In the present effective 
theory this aspect is reflected by two amplitude degrees 
of freedom present in (Q). More detailed discussion of 
these issues is given in Refs. [ ^Ops] ]. 

We have thus arrived at an effective theory of a spin- 
charge separated system containing one phase degree of 
freedom f2 and two amplitudes, pa and pb- Deep in the 
superconducting phase, where both amplitudes arc finite, 
the physics of (Q) will be very similar to that of a conven- 
tional GL theory. In the situations where the supercon- 
ducting order parameter is strongly suppressed, such 
as in the vortex core, near an impurity or a wall, the new 
theory has an extra degree of richness, associated with 
the fact that it is sufficient (and generally preferred by 
the energetics) when only one of the two amplitudes is 
suppressed. Since the two amplitudes play very different 
roles in the electronic excitation spectrum, the effective 
theory (0) will lead to a number of nontrivial effects. 

To illustrate this consider what will happen in the core 
of a superconducting vortex. Under the influence of the 
magnetic field the phase £1 will develop a singularity such 
that V51 ~ 1/r close to the vortex center. For the free 
energy to remain finite the amplitude prefactor in the 
second term of Eq. (Hj must vanish for r — > 0. This is 
analogous to | M/ 1 vanishing in the core of a conventional 
vortex. In the present case, however, it is sufficient when 
the product p^Pb vanishes. Since suppressing any of the 
two amplitudes costs condensation energy, in general only 
one amplitude will be driven to zero. Which of the two 
is suppressed will be determined by the energetics of the 
amplitude term (^). On general grounds we expect that 
the state in the vortex core will be the same as the cor- 
responding bulk "normal" state obtained by raising tem- 
perature above T c . Thus, very crudely, we expect that 
holon vortex will be stable in the underdoped while the 
spinon vortex will be stable in the overdoped region of 
the phase diagram Figure |l|. 

An important point by which our approach differs from 
the SNL theory is that in the present theory both types 
of vortices carry the same superconducting flux quantum 
$o and thus compete on equal footing. This is a direct 
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consequence of our assumption of the vanishing phase 
stiffness a. 

In what follows we study in detail the vortex solutions 
of the free energy (^). Our main objective is to obtain 
the precise estimates for the energy of the two types of 
vortices as a function of temperature and doping and de- 
duce the corresponding phase diagram for the state inside 
the vortex core. We show that for generic parameters in 
(Q) the singly quantized holon vortex with a pseudogap 
spectrum in the core can be stabilized over a large por- 
tion of the superconducting phase, as required by the 
experimental constraints discussed above. 



II. SOLUTION FOR A SINGLE VORTEX 

A. General considerations 

In order to provide a more quantitative discussion we 
now adopt some assumptions about the coefficients en- 
tering the free energy ([!]). We assume that 



Ti = oti(T - T{), 



6, A, 



(7) 



where 



Pi 



Ps 



4pi 



(12) 



is the effective superfluid density. The term in brack- 
ets can be identified as twice the conventional superfluid 
velocity 

v s = -X7Q, — eA. 
2 

Making use of the Ampere's law 47rj = V x B we see that 
Eq. (fTy) specifies the supercurrent in terms superfluid 
density and velocity: j = 2ep s v s . Minimization of (Q) 
with respect to then implies V • j = 0; the supercurrent 
is conserved. 

Minimizing the free energy (Q) with respect to the am- 
plitudes results in the pair of coupled GL equations: 



V 2 



PA + TAP A 



UAP\ + VplpA 



^PapI 



Pi) 2 



o, 



(13a) 



where T; are corresponding "bare" critical temperatures, 
which we assume depend on doping concentration x in 
the following way: 



Ta = T Q (2x r , 



T b = T x. 



(8) 



Here x m denotes the optimal doping and Tq sets the over- 
all temperature scale. We furthermore assume that ut 
and v are all positive and independent of doping and 
temperature. It is easy to see that such choice of param- 
eters qualitatively reproduces the bulk phase diagram of 
cuprates in the x-T plane shown in Figure [3]. The effect 
of the w-term is to suppress T c from its bare value away 
from the optimal doping. In real systems fluctuations 
will lead to additional suppression of T c which we do not 
consider here. 

In the absence of perturbations the bulk values of the 
amplitudes are given by 



Pa = ~(rAUb - r b v)/D, 
-(nuA - tav)/D, 



Pi 



(9) 



with D = UfcMA — v 2 . In analogy with conventional GL 
theories we may define coherence lengths for the two am- 
plitudes Go) 



Ca 2 = _ ( r A - nv/u b ), 

£b~ 2 = ~(n - tav/ua), 



(10) 



one of which always diverges at T c as (T — Tc)^ 1 / 2 . 

Minimization of the free energy (Q) with respect to the 
vector potential A yields an equation 



V x V x A = ep s (Vtt - 2eA), 



(11) 



V 2 pb + r b p b + u b pl + vp\p b + 



(4pi + P 2 ) 



2\2 s 



(13b) 



We are interested in the behavior of the amplitudes in 
the vicinity of the vortex center. In this region, for a 
strongly type-II superconductor, we may neglect the vec- 
tor potential A in the superfluid velocity v s . In a singly 
quantized vortex Q winds by 2tt around the origin lead- 
ing to a singularity of the form v s ~ i Vfi = ip/2r. First, 
for the holon vortex we assume that p b vanishes in the 
core as some power p b (r) ~ r v and Pa(j") ~ Pa remains 
approximately constant. Eq. (13t) then becomes 



/I 2\ u-2 

(4 ~ v > r 



(r b + vp 2 A )r" 



0, (14) 



where we have neglected p\(r) compar ed t o 4p A in the 
denominator of the last term in Eq. (13t). The most 
singular term in Eq. (|lj) is the first one and we must de- 
mand that the coefficient of r v ~ 2 vanishes. This implies 
v = i. The asymptotic short distance behavior of the 
holon amplitude therefore can be written as 



Pb(r) ~ c b p b 



6 



1/2 



(15) 



where c b is a constant of order unity which may be de- 
termined by the full integration of Eqs. (|i"3|). Similar 
analysis of Eq. (13a) in the vicinity of the spinon vortex 
yields 



p A (r) ~ cap a 



(16) 
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with pb approximately constant. 

We notice the different power laws in the holon and 
spinon results. Operationally this difference arises from 
different numerical prefactors of the respective superfluid 
velocity terms in Eqs. (^). Physically, the unusual r 
dependence of the holon amplitude in the core reflects 
the fact that the field b describes a condensate of single 
holons, each carrying charge e. Superconducting vortex 
with the flux quantum <&o represents a magnetic "half- 
flux" for the holon field which results in non-analytic be- 
havior of Pb(r) at the origin. Singly quantized holon vor- 
tex is therefore a peculiar object and we shall discuss it 
more fully in Section III. Here we note that the physical 
superconducting order parameter amplitude \^\ = pap\ 
remains analytic in the core of both the spinon and the 
holon vortex. 




B. Holon vs. spinon vortex: the phase diagram 

We are now in the position to estimate the energies 
of the two types of vortices and deduce the phase dia- 
gram for the "normal" state in the vortex core. To this 
end we consider a single isolated vortex centered at the 
origin. The total vortex line energy can be divided into 
electromagnetic and core contributions J3^|. The elec- 
tromagnetic contribution consists of the energy of the 
supercurrents and the magnetic field outside the core re- 
gion. It may be estimated by assuming that the ampli- 
tudes pa and pb have reached their bulk values Pa and 
pb respectively. Taking curl of Eq. (^) and noting that 
V x V£l = 27r<5(r) for a singly quantized vortex we obtain 
the London equation for the magnetic field B = V x A 
of the form 



FIG. 2. Vortex core phase diagram for GL parameters cho- 
sen as follows: a A = 0.13, a b = 0.10, T = 200K, x m = 0.2, 
ua = U), — 1.0 and v = 0.5. Dashed line marks the phase 
boundary T g {x) obtained from Eq. ( pB] ) while the solid cir- 
cles correspond to the numerical calculation with the same 
parameters. 



while the other one stays constant and equal to its bulk 
value. This is a very crude approximation which we jus- 
tify below by an exact numerical computation. With 
these assumptions, the core energy is 



E (i) 



( $o 



(21) 



where i — A, b for spinon and holon vortex respectively 
and 



B - \ 2 V 2 B = <f> 5(r) 



where 



A = 87re z 



4pl 



Pi 



(17) 



(18) 



has the meaning of the London penetration depth for 
the effective GL theory (|J) . Aside from the unusual form 
of A, Eq. (|l7|) is identical to the conventional London 
equation. The corresponding electromagnetic energy is 
therefore the same for both types of vortices and can be 
calculated in the usual manner |52 33 3(]] obtaining 



^ em "(5a-) inK ' 



(19) 



with k = A/ max(£A, £f>) being the generalized GL ratio. 

To estimate the core contribution to the vortex line en- 
ergy we assume that one of the amplitudes is suppressed 
to zero in the core 



Pi(r) = 0, r < & 



(20) 



A; 



(22) 



Such a crude approximation overestimates the core en- 
ergy. A more accurate analysis j$2 33|, which we do not 
pursue here, allows for a more realistic variation of pi (r) 

in the core and indicates that the value of -Ecore has the 
same form as Eq. ( ^l| ) multiplied by a numerical factor 
Ci m 0.5 |3{]j37]]. Thus, the total energy of the vortex line 
can be written as 



E (i) 



(—)' 
\4ttXJ 



In k + Ci 



V47rAi 



(23) 



where again i = A, b for spinon and holon vortex respec- 
tively. Eq. (H|) parallels the Abrikosov expression for the 
vortex line energy in a conventional GL theory ]32] ] where 
A and Xi are identical and equal to the ordinary London 
penetration depth. 

In the vortex state described by the free energy (Q) 
the vortex with lower energy E^ 1 ' will be stabilized. Eq. 
( p3| ) implies that the difference in energy between the two 
types of vortices comes primarily from the core contribu- 
tion, as expected on the basis of the physical argument 
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presented above. Condition Aa = Ab marks the tran- 
sition point between the two solutions. For fixed GL 
parameters Tq, x m , ctj, it, and v this defines a transition 
line in the x-T plane. According to (^2| ) the equation for 
this line is 



Pa{x,T) = p b (x,T). 



(24) 



Using Eqs. (|7]-|9|) one can obtain an explicit expression 
for the transition temperature T g between two types of 
vortices as a function of doping 



with 



T g (x) = T 



1-/5 \-l3~ 1 



Otb(uA 



(25) 



(26) 



Eq. (|25| ) describes a straight line in the x-T plane, orig- 
inating at [x m ,Tox m ], i.e. maximal T c at optimum dop- 
ing, and terminating at [2x m /(l + /3),0]. Gcnerically, we 
expect that parameters on and Ui will be comparable in 
magnitude for the holon and spinon channels. Parameter 
j3 defined in Eq. ([26j) will therefore be of order unity. The 
typical situation for (3 = 0.77 is illustrated in Figure |^. 
More generally the quartic coefficients Ui and v could ex- 
hibit weak doping and temperature dependences leading 
to a curvature in the phase boundary. 

The appealing feature of the present theory is that pa- 
rameter (3 may vary from compound to compound. Thus, 
the experimental fact that in BSCCO the pseudogap in 
the core persists into the overdoped region is easily ac- 
counted for in the present theory. It would be interesting 
to see if the transition from holon to spinon vortex as 
a function of doping could be experimentally observed. 
A good candidate for such observation would be LSCO, 
where the transport measurements in pulsed magnetic 
fields Q| established a metal-insulator transition around 
optimal doping, i.e. (3 ~ 1. The current theory predicts 
a holon vortex with the pseudogap spectrum in the un- 
derdoped (insulating) region and spinon vortex with con- 
ventional metallic spectrum on the overdoped side. 



C. Numerical results 

In order to put the above analytical estimates on firmer 
ground we now pursue numerical computation of the vor- 
tex line energy. For simplicity we consider the strongly 
type-II situation (n 3> 1) where the vector potential term 
in v s can be neglected to an excellent approximation, as 
long as we focus on the behavior close to the core. We are 
then faced with the task of numerically minimizing the 
free energy (|j) with respect to the two cylindrically sym- 
metric amplitudes pa{t) aud Pb(f)- As noted by Sachdev 
pof direct numerical minimization of the free energy (EJ) 



(a) holon vortex 
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FIG. 3. Order parameter amplitudes near a single isolated 
vortex for GL parameters specified in Figure |^. The holon 
vortex is plotted for T — and x — 0.22 (implying coherence 
lengths £a = 0.63 and £b = 0.70), while the spinon vortex 
is plotted for T = and x = 0.24 (implying £a = 0.75 and 
£b = 0.60). 



provides a more robust solution than the numerical inte- 
gration of the coupled differential equations (|13| ). 

We discretize the free energy functional M) on a disk 
of a radius R 3> £j in the radial coordinate r with up 
to N = 2000 spatial points. We then employ the Polak- 
Ribiere variant of the Conjugate Gradient Method |3§| ] 
to minimize this discretized functional with respect to 
PA(fj) and Pb(fj), initialized to suitable single vortex 
trial functions. The procedure converges very rapidly 
and the results are insensitive to the detailed shape of 
the trial functions as long as they saturate to the correct 
bulk values outside the vortex core. 

Typical results of our numerical computations are dis- 
played in Figure (||) and are in complete agreement with 
the analytical considerations of the preceding subsec- 
tions. Note in particular that Pb{f) in the holon vor- 
tex vanishes with infinite slope, consistent with Eq. (|l5|). 
Plotting p\{r) confirms that the exponent is indeed 1/2. 
In the spinon vortex Pa(?) is seen to vanish linearly as 
expected on the basis of Eq. (pit). The nonvanishing or- 
der parameter is slightly elevated in the core reflecting 
the effective "repulsion" between the two amplitudes con- 
tained in the v-term of the free energy. The results for 
the spinon vortex are consistent with those of Ref. |3(J] . 

We explored a number of other parameter configura- 
tions and obtained similar results. We find that Eq. (BJ) 



G 



is a good predictor of the transition line between the 
holon and spinon vortex, although the precise numeri- 
cal value of the transition temperature T g for given x 
tends to deviate slightly from the value predicted by Eq. 
(p5|). This is illustrated in Figure (^J) where we compare 
the vortex core phase diagrams obtained numerically and 
from Eq. (p5|). Interestingly, the deviation always tends 
to enlarge the holon vortex sector of the phase diagram 
at the expense of the spinon vortex sector. This is pre- 
sumably because the sharper ~ y/r suppression of the 
holon order parameter in the core costs less condensa- 
tion energy. 



/gl = /amp + pi(V0 - 2a) 2 + p 2 b (V9 - a - eA) 2 
+ i-(Vx A) 2 + |(Vxa) 2 , (27) 

with /amp specified by Eq. (|J). Now consider situation 
in which the sample is subjected to uniform magnetic 
field B = V x A. Two scenarios (discussed previously 
by SNL) appear possible. In the first, the internal gauge 
field develops no net flux, (V x a) = 0, and the holon 
phase 9 develops singularities in response to A such that 



V x V0 = 27r^<5(r 



III. GAUGE FLUCTUATIONS AND THE 
SPECTRAL PROPERTIES IN THE CORE 

Theory of the vortex core based on the effective ac- 
tion (|j) appears to yield results consistent with the STM 
data on cuprates in that it implies stable holon vor- 
tex solution over the large portion of the superconduct- 
ing phase diagram. The state inside the core of such a 
holon vortex is characterized by vanishing amplitude of 
the holon condensate field, |6| =0, and a finite spin gap 
|A| m Abuik- This is the same state as in the pseudogap 
region above T c . One would thus expect the electronic 
spectrum in the core to be similar to that found in the 
normal state of the underdoped cuprates, in agreement 
with the data ||J|. The holon vortex with this prop- 
erty carries conventional superconducting flux quantum 
$o, in accord with experiment. This general agreement 
between theory and experiment would suggest that the 
effective action (|J) provides the sought for phenomeno- 
logical description of the vortex core physics in cuprates. 
In what follows we amplify our argumentation that it is 
also tenable in a broader theoretical context in that it 
naturally follows from the U(l) slave boson models ex- 
tensively studied in the classic and more recent high-T c 
literature. We then provide a more detailed discussion of 
the vortex core spectra and propose an explanation for 
the experimentally observed core bound states. 



A. Significance of the / gaU gc term 

Derivation of the effective action 0) from the more 
general U(l) action (^|) hinges on our assumption that the 
stiffness a of the gauge field a is low and that the /gauge 
term (||) can be neglected. Assumption of large <j by 
SNL leads to very different vortex solutions |30],[3l) which 
appear inconsistent with the recent experimental data. 
We first expand on our discussion as to why is /gauge term 
important and then we argue why it may be permissible 
to neglect it in the realistic models of cuprates. 

To facilitate the discussion let us rewrite Eq. (||) by 
resolving the complex matter fields into amplitude and 
phase components: 



where rj denotes the vortex positions. The holon ampli- 
tude pb is driven to zero at Vj , essentially to prevent the 
free energy from diverging due to the singularity in the 
phase gradient. Since holons carry charge e, each vor- 
tex is threaded by flux hc/e, i.e. twice the superconduct- 
ing flux quantum $o = hc/2e. This solution represents 
the doubly quantized holon vortex lattice, considered by 
SNL. 

In the second scenario a develops a net flux such that 
a ps — eA, which screens out the A field in the holon term 
but produces a net flux — 2eA in the spinon term. In re- 
sponse to this flux, spinon phase <f> develops singularities 
such that 



V x = 2ttJ2 S ( 1 



corresponding to the spinon vortex lattice. Tj denotes 
vortex positions which will be different from Tj since at 
the fixed field B there will be twice as many spinon vor- 
tices as holon vortices. (Spinon vortices carry conven- 
tional superconducting quantum of flux $o-) I n this case 
Pa is driven to zero at the vortex centers. In this scenario 
one pays a penalty for nucleating the net flux in V x a 
due to last term in Eq. (|27]). This energy cost can be 
estimated as 



\4ttX J 



(28) 



per vortex. Stiffness a must be small enough so that E a 
is small compared to the vortex energy (E3). Taking the 
dominant Eem term and neglecting In k this implies that 



cr < 



f 



8ttc 2 



(29) 



which is the same condition as considered in Ref. [pPf . 

Now consider a third scenario in which a singly quan- 
tized holon vortex emerges. As a starting point consider 
the spinon vortex solution just described. In the un- 
derdoped regime the amplitude piece / am p would favor 
suppressing the holon amplitude in the core instead of 
the spinon amplitude but according to our previous con- 
siderations this would ordinarily require formation of a 
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doubly quantized vortex whose magnetic energy is too 
large. However, if the gauge field stiffness a is sufficiently 
small, the system could lower its free energy by setting 
up singularities in a which would precisely cancel the sin- 
gularities in V</> and shift them to the holon term. To 
arrive at this situation imagine contracting the initially 
uniform flux V x a so that it becomes localized in the 
individual vortex core regions. Taking this procedure to 
the extreme, i.e. taking the limit a — > 0, the gauge field 
will form "flux spikes" of the form 



2(V x a) = -V x = -2tt ^ <5(i 



(30) 



completely localized at the vortex centers. Gauge field 
of this form indeed completely cancels the singularities 
in the spinon phase gradient in Eq. ( |2^ ) and pa is no 
longer forced to vanish in the core. The singularities now 
appear in the holon term, but they stem from a rather 
that V# which remains nonsingular. Consequently, pb is 
forced to vanish in the vortex cores. By construction the 
vortices are located at ij and are therefore singly quan- 
tized. This is the singly quantized holon vortex discussed 
in the framework of the free energy ([|). Based on the 
above discussion the singly quantized holon vortex can 
be thought of as a composite object formed by attaching 
half quantum (h/2) of the fictitious gauge flux V x a to 
the spinon vortex. Within the full compact U(l) theory 
this is essentially equivalent to the Z2 vortex discussed by 
Wen J|9| in the framework of topological orders in spin 
liquids. 

In the framework of the free energy ( p7| ) one pays a 
penalty for such a singular solution due to the gauge stiff- 
ness term. In the present continuum model this penalty 
per single vortex is actually infinite, since according to 
Eq. ( |3(i| ) it involves a spatial integral over [5(r — ij)] 2 . 
Thus, in the continuum model the singular solutions of 
this type are prohibited. In reality, however, we have to 
recall that our effective action (||) descended from a mi- 
croscopic lattice model for spinons and holons in which 
the gauge field a lives on the nearest neighbor bonds of 
the ionic lattice. The ionic lattice constant d therefore 
provides a natural short distance cutoff and the delta 
function in Eq. (^) should be interpreted as a flux quan- 
tum $0 piercing an elementary plaquette of the lattice. 
The energy cost per vortex thus becomes finite and is 
given by 



ae 2 f $ 

E °- — It 



(31) 



Again, for the solution to be stable, E'„ must be negligible 
compared to the vortex energy (|23j). This implies 



f 



8ir 2 e 2 



(32) 



which is a much more stringent condition than (|29|) since 
in cuprates d<A. 



When condition (|32|) is satisfied it is permissible to ne- 
glect the /gauge term in the effective action (|^) and it 
becomes fully equivalent t o (Eh as far as the vortex solu- 
tions are concerned. Eq. (p2|) gives the precise meaning 
to the requirement of the weak stiffness of the gauge field 
loosely stated when deriving the effective action (|j) . 



B. Microscopic considerations 

As mentioned in the introduction, the gauge field a 
has no dynamics in the original U(l) microscopic model, 
as it only serves to enforce a constraint on spinons and 
holons. The stiffness term (|J) in the effective theory was 
assumed to arise in the process of integrating out the mi- 
croscopic degrees of freedom |3(],[jl| . While such term is 
certainly permitted by symmetry, assessing its strength a 
is a nontrivial issue since even deep in the superconduct- 
ing phase neither holons nor spinons are truly gapped. 
Thus, in general, integrating out these degrees of freedom 
may lead to singular and nonlocal interactions between 
the condensate and the gauge fields. To our knowledge 
the procedure has not been explicitly performed for the 
U(l) model and the precise form or magnitude of the 
gauge stiffness term is unknown. General considerations 
|35| suggest that the gauge stiffness term is negligible in 
the class of models with exact local U(l) symmetry con- 
necting the phases of holons and spinons. 

Consider now an intermediate representation of the 
problem where only high energy microscopic degrees of 
freedom have been integrated out. In the presence of 
a cutoff this is a well defined procedure even for gap- 
less excitations, as explicitly shown by Kwon and Dorsey 
pof for a simple BCS model. The corresponding effec- 
tive Lagrangian density of the present U(l) model can be 
written as 



C, = if(^0-2a Al ) 2 + ^(9 / 



^■^-fx) /amp 



2 p/ 2 

+ {d^ - 2o„) J& + (d„e - a„ - eA„) J£ 

+ C sp [tp sp ,ipl p ;p A } + Ch[i>hA h \ Pb] +£em[Au]- ( 33 ) 

The Greek index p runs over time and two spatial dimen- 
sions, re? are compressibilities of the holon and spinon 
condensates, while 



A, 6, 



J 



1,2, 



(34) 



are the respective phase stiffnesses. Jg p and j£ are spinon 
and holon three currents respectively and C sp and Ch are 
the low energy effective Lagrangians for the fermionic 
spinon field ip sp and bosonic holon field ip^. £em is 
the Maxwell Lagrangian for the physical electromagnetic 
field. Thus, £ e ff describes an effective low energy theory 
of spinons and holons coupled to their respective collec- 
tive modes and a fluctuating U(l) gauge field. Similar 
theory has been recently considered by Lee El . 
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The precise form of the microscopic Lagrangians £ sp 
and Lh is not important for our discussion. The salient 
feature which we exploit here is that only the amplitude 
of the respective condensate field enters into C sp and £ n . 
Coupling to the phases and the gauge field is contained 
entirely in the Dopplcr shift terms [second line of Eq. 
(|33|)]. Such form of the coupling is largely dictated by the 
requirements of the gauge invariance and the particular 
form Eq. (33) can be explicitly derived by gau ging a way 
the respective phase factors from the ip fields [£9|,[l0| . 

The gauge field a M enters the effective Lagrangian ( pjj ) 
only via two gauge invariant terms: (d^cf) — 2a^) and 



- eAfj), which may be interpreted as the three 



velocities of the spinon and holon condensates respec- 
tively. Furthermore, the only coupling between holons 
and spinons arises from a M . Therefore, if we now pro- 
ceed to integrate out the remaining microscopic degrees 
of freedom from £ e ff, the two velocity terms will not mix. 
This consideration suggests that upon integrating out all 
of the microscopic degrees of freedom, the resulting gauge 
stiffness term will be of the form 



f 

■I da 



2 



V x (2a- V(/>)] 2 



?[Vx(a- 



eA- V6>)] 5 



(35) 



Clearly, such term is permitted by the gauge symme- 
try. Furthermore, we note that for smooth (i.e. vortex 
free) configurations of phases the gradient terms will con- 
tribute nothing and we recover the gauge term considered 
in Ref. Q. 

In the presence of a vortex in 4> or 9 the fg augc term 
will contribute formally divergent energy. Regularizing 
this on the lattice, as discussed above Eq. ([H]), this en- 
ergy will become finite and can be interpreted simply 
as the energy of the spinon or holon vortex core states, 
which have been integrated out. In the microscopic the- 
ory ( |33| ) such energy would arise upon solving the rele- 
vant fermionic or bosonic vortex problem. 

We stress that, as concluded in the preceding subsec- 
tion, the main theoretical obstacle to the formation of a 
singly quantized holon vortex in the original SNL theory 
was the appearance of a formally divergent contribution 
in the / gaU gc term (H). The argument above sug gest s 
that /gauge in Eq. ([|) should be replaced by Eq. (35), 
in which such formally divergent contribution appears 
for arbitrary vortex configuration and upon regulariza- 
tion has a simple physical interpretation in terms of the 
energy of the vortex core states. Usage of the physically 
motivated term (35) in place of (||) therefore removes the 
bias against the singly quantized holon vortex solution, 
which appears to be realized in real materials. With ( |35| ) 
any bias between the holon and spinon vortex solutions 
can result only from the difference between the two stiff- 
ness constants a a and ov It is reasonable on physical 
grounds to assume that constants <ta and Ub are of the 
similar magnitudes. Furthermore, on the basis of Ref. 
p5| we expect these constants to be negligibly small in 



the physically relevant models. Consequently we expect 
that neglecting the /gauge term as in our derivation of 
effective action ([!]) will result in accurate determination 
of the phase diagram for the state in the vortex core. 



C. Vortex core states 

The phenomenological theory based on the effective 
action (|J) does not allow us to address the interesting 
question of the nature of the fermionic states in the vor- 
tex core. To do this we need to consider the microscopic 
Lagrangian density (|3^). While the fully self consistent 
calculation is likely to be prohibitively difficult, one can 
obtain qualitative insights by first solving the GL the- 
ory %) as described in Sec. II, and then using the order 
parameters /?a and pb as an input to the fermionic and 
bosonic sectors of the theory specified by Eq. (|3^) . The 
work on a detailed solution of this type is in progress. 
Here we wish to point out some interesting features of 
such a theory and argue that it may indeed exhibit struc- 
ture in the low energy spectral density similar to that 
found experimentally ||||. 

It is instructive to integrate out the gauge fluctuations 
from the Lagrangian ( |3^ ) as first discussed by Lee [ p0[ . 
Since £ e ff is quadratic in the integration can be ex- 
plicitly performed resulting in the Lagrangian of the form 

C' eB = ^M<) 2 - /amp + £ E M 

2K b I n jfA \ , i^A i n 

<m A -t- K b 



4< + k£ 



l l 



2 + Kt 



(2J& + (36) 



where = 4k^k^/(4k^ + and 



v» = {d,9--d^-eA,) 



(37) 



is the physical superfluid velocity. The first line re- 
produces the GL effective action (^) for the condensate 
fields, the second line describes the Doppler shift cou- 
pling of the superfluid velocity to the microscopic cur- 
rents, and the third line contains spinon and holon pieces 
with additional current-current interactions generated by 
the gauge fluctuations (2(J . 

We now discuss the physical implications of Eq. (|3^ ) 
for the two types of vortices. We focus on the static 
solutions (i.e. we ignore the time dependences of various 
quantities, e.g. taking v° s = 0) of C' eS in the presence of 
a single isolated vortex. We are interested in the local 
spectral function of a physical electron. This is given 
by a convolution in the energy variable of the spinon 
and holon spectral functions. According to the analysis 
presented in Ref. |2q] , at low temperatures the electron 
spectral function will be essentially equal to the spinon 
spectral function. Convolution with the holon spectral 







function which is dominated by the sharp coherent peak 
due to the condensate merely leads to a small broadening 
of the order T. In the following we therefore focus on the 
behavior of spinons in the vicinity of the two types of 
vortices. 

By inspecting Eq. (|36| ) it is easy to see that the ex- 
citations inside the spinon vortex will be qualitatively 
very similar to those found in the conventional vortex de- 
scribed by the weak coupling <i-wave BCS theory [p|-^0[. 
In particular according to Eq. (|l6|) we have ka ~ r, 
and Kb ~ const in the core. Recalling furthermore that 
|v s | ~ 1/rwe observe that the spinon current J sp is cou- 
pled to a term that diverges as 1/r in the core (just as in 
a conventional vortex), while the holon current J/j is cou- 
pled to a nonsingular term. Thus, one may conclude that 
holons remain essentially unperturbed by the phase sin- 
gularity in the spinon vortex while the spinons obey the 
essentially conventional Bogoliubov-de Gennes equations 
for a d-wave vortex. 

In the holon vortex the situation is quite different. Ac- 
cording to Eq. ( [l5| ) we have Kb ~ r and k\ ~ const in the 
core. The spinon current J sp is now coupled to a nonsin- 
gular term (1/r divergence in v s is canceled by Kb ~ r). 
Therefore, there will be no topological perturbation in 
the spinon sector and we expect the spinon wavefunctions 
to be essentially unperturbed by the diverging superfluid 
velocity. Spinon spectral density in the core should be 
qualitatively similar to that far outside the core. This is 
our basis for expecting a pseudogap-like spectrum in the 
core of a holon vortex. 

We now address the possible origin of the experimen- 
tally observed vortex core states within the present 
scenario for a holon vortex. To this end consider the 
effect of the last term in Eq. ( |36| ) which we ignored so 
far. Upon expanding the binomial the temporal compo- 
nent is seen to contain a density-density interaction of 
the form </ s p J/j where is the local density of uncon- 
densed holons. Since the holon order parameter vanishes 
in the core and the electric neutrality dictates that the 
total density of holons must be approximately constant 
in space, we expect that uncondesed holon density will 
behave roughly as 

J°(r) = Pb~ Pb{r); 

J°(r) will have a spike in the core of a holon vortex. Inso- 
far as (r) can be viewed as a static potential acting on 
spinons, the uncondensed holons in the vortex core can be 
thought of as creating a scattering potential, akin to an 
impurity embedded in a rf-wave superconductor. In fact, 
formally the spinon problem is identical to the problem 
of a fermionic quasiparticle in a <i-wave superconductor 
in zero field in the presence of a localized impurity po- 
tential. It is known that such problem exhibits a pair 
of marginally bound impurity states [Uj at low energies 
which result in sharp resonances in the spectral density 
inside the gap. Such states have been extensively stud- 
ied theoretically [}42|-[45fl and their existence was recently 



confirmed experimentally by Pan et al. pL We propose 
here that, within the formalism of Eq. (36|), the same 
mechanism could give rise to the low energy quasiparti- 
cle states in the core of a holon vortex. Such structure, 
if indeed confirmed by a microscopic calculation, could 
explain the spectral features observed experimentally in 
the vortex cores of cuprate superconductors MM . 



IV. CONCLUSIONS 

Scanning tunneling spectroscopy of the vortex cores 
affords a unique opportunity for probing the underly- 
ing "normal" ground state in cuprate superconductors. 
The existing experimental data on YBCO and BSCCO 
strongly suggest that conventional mean field weak cou- 
pling theories f7]-|l3| fail to describe the physics of the 
vortex core. Our main objective was to develop a the- 
oretical framework for understanding these spectra and 
the nature of the strongly correlated electronic system 
which emerges once the superconducting order is sup- 
pressed. We have shown that phenomenological model 
(||) based on a variant of the U(l) gauge field slave bo- 
son theory [M contains the right physics, provided that 
the gauge field stiffness is vanishingly small. The latter 
assumption is consistent with the general arguments in- 
volving local gauge symmetry p5j . In such a theory the 
gauge field can be explicitly integrated out, resulting in 
the effective action (^) which contains one phase degree 
of freedom representing the phase of a Cooper pair and 
two amplitude degrees of freedom representing the holon 
and spinon condensates. 

Analysis of the effective theory (|J) in the presence of a 
magnetic field establishes existence of two types of vor- 
tices, spinon and holon, with contrasting spectral proper- 
ties in their core regions. Our holon vortex is singly quan- 
tized and therefore differs in a profound way from the 
doubly quantized holon vortex discussed by SNL |5(i|,[n| . 
As indicated in Figure |^ such a singly quantized holon 
vortex is expected to be stable over the large portion of 
the phase diagram on the underdoped side. Quasiparti- 
cle spectrum in the core of a holon vortex is predicted to 
exhibit a "pseudogap" , similar to that found in the un- 
derdoped normal region above T c . This is consistent with 
the data of Renner et al. who pointed out a remark- 
able similarity between the vortex core and the normal 
state spectra in BSCCO. Spinon vortex, on the other 
hand, should be virtually indistinguishable from the con- 
ventional d-w&ve BCS vortex and is expected to occur on 
the overdoped side of the phase diagram. Transition from 
the insulating holon vortex to the metallic spinon vortex 
as a function of doping is a concrete testable prediction 
of the present theory. 

Phenomenological theory based on the effective action 
([|) does not permit explicit evaluation of the electronic 
spectral function. To this end we have considered the cor- 
responding microscopic theory (pq) and concluded that 
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holon vortex will indeed exhibit a pseudogap like spec- 
trum. Such qualitative analysis furthermore suggests a 
plausible mechanism for the sharp vortex core states ob- 
served in YBCO § and BSCCO §. We stress that con- 
ventional mean field weak coupling theories yield neither 
pseudogap nor the core states. In the core of a holon vor- 
tex such states will arise as a result of spinons scattering 
off of the locally uncondensed holons, in a manner anal- 
ogous to the quasiparticle resonant states in the vicinity 
of an impurity in a e?-wave superconductor [^l]-^5| . The 
latter conclusion is somewhat speculative and must be 
confirmed by explicitly solving the fermionic sector of 
the microscopic theory (|36|). 

On a broader theoretical front the importance of the 
vortex core spectroscopy as a window to the normal state 
in the T — > limit lies in its potential to discriminate 
between various microscopic theories of cuprates. It is 
reasonable to assume that the observed pseudogap in the 
vortex core reflects the same physics as the pseudogap 
observed in the normal state. This means that the mech- 
anism responsible for the pseudogap must be operative 
on extremely short lengthscales, of order of several lattice 
spacings. The U(l) slave boson theory considered in this 
work apparently satisfies this requirement. Obtaining the 
correct vortex core spectral functions could serve as an 
interesting test for other theoretical approaches describ- 
ing the physics of the underdoped cuprates J28|,[l6],[l7| . 

It will be of interest to explore the implications of the 
effective theories ([!]) and ( |36| ) in other physical situa- 
tions. Of special interest are situations where the holon 
condensate amplitude is suppressed, locally or globally, 
giving rise to "normal" transport properties (vanishing 
superfluid density) but quasiparticle excitations that are 
characteristic of a superconducting state. These include 
the spectra in the vicinity of an impurity, twin boundary 
or a sample edge. In the latter case one might hope to ob- 
serve a signature of the zero bias tunneling peak anomaly 
(normally seen for certain geometries deep in the super- 
conducting phase in the optimally doped cuprates) even 
above T c in the underdoped samples. 
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Note added in proof. After submission of this 
manuscript we learned about complementary microscopic 
treatments of the spin-charge separated state in the vor- 
tex core within U(l) @ and SU(2) |§ slave boson 
theories. The former agrees qualitatively with our phe- 
nomenological theory. Ref. |49| proposes a new type 



of vortex which takes advantage of the larger symme- 
try group SU(2). In a related development Senthil and 
Fisher discussed a Z2 vortex (which is essentially 
equivalent to our singly quantized holon vortex) and pro- 
posed a "vison detection" experiment based on trapping 
such a vortex in the hole fabricated in a strongly under- 
doped superconductor. Here we wish to point out that 
the experiment will produce the same general outcome in 
a system described by the U(l) theory where the role of 
a vison will be played by a flux quantum of the fictitious 
gauge field a. 
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